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ABSTRACT 


The natural convection heat transfer from an array of heaters flush mounted on a 
vertical test surface in response to both step and periodic input power has been investi- 
gated. Baseline surface temperature measurements in water were obtained for several 
steady power inputs, ranging from 0.2 W to 1.0 W for the vertical test surface. Also, the 
effects of the presence of a parallel shroud near the test surface on measured temperature 
values were studied. The effects of periodic variations in input power on the heater ele- 
ment temperatures were then examined at selected locations. The resulting transient 
heater temperatures were measured for various amplitudes, frequencies and average lev- 
els of the periodic input power. These transient measurements are compared with a 
two-dimensional numerical simulation of the relevant fluid flow and heat transfer proc- 
esses. The computed heater surface temperatures supported the trends of the measure- 
ments and were within the expected range. 
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I. INTRODUCTION 


A considerable increase in the volumetric power generation rates in clectronic com- 
ponents has resulted due to both continued growth in functional capability and rapid 
decrease in feature size. The future enhancements in speed and size reduction of elec- 
tronic systems are constrained by available heat removal techniques. 

The performance and reliability of microelectronic components are Very sensitive to 
junction temperature. The average failure rates of chips increase exponentially as the 
junction temperature increases Bar-Cohen and kraus [Ref. 1]. In recent years, due to 
the relatively high heat transfer rates and other potential advantages such as high reli- 
ability, low cost, reduced noise and simpler design, complete immersion of microelec- 
tronic circuits in dielectric liquids has been examined. 

Relatively few studies on natural convection liquid immersion cooling are available 
in literature. The heat transfer characteristics of flush mounted heat sources on vertical 
surface in several liquids were studied by Baker et al [Ref. 2]. Due to the three dimen- 
sional effects, large increases in heat transfer were noticed. Both flush and protruding 
arrangements were investigated in water and R-113 by Park and Bergles (Ref. 3]. Cor- 
relations for natural convection in a vertical channel dissipating uniform heat flux from 
sides were presented dy Bar-Cohen and Schweitzer [Ref. 4]. Flow visualization and heat 
transfer measurements for natural convection from a long heated protrusion on a verti- 
cal insulated wall of a rectangular chamber filled with water were carried out by Kelleher 
and Knock [Ref. 5]. These measurements were supported by a computational investi- 


gation by Lee et al [Ref. 6]. 








Flow visualizations and component surface temperature measurements for natural 
convection cooling of a 3 by 3 array of discrete heated frotrusions on the vertical wall 
of a rectangular enclosure filled with a dielectric fluid were performed by Joshi et al 
[Ref. 7]. In other experimental studies, the heat transfer and flow characteristics of 
protruding heat sources on a vertical surface (Ref. 8] and vertical channel [Ref. 9] were 
studied by Joshi et al. Flow visualizations and component surface temperature meas- 
urements Were carried out for several power dissipation levels and channel spacings. For 
the smallest channel spacing, a significant increase in the component surface temper- 
ature due to a reduction in fluid velocity was observed. 

Natural convection adjacent to a three column array of fifteen per column uniformly 
heated components flush mounted on a vertical flat plate in water was examined by 
Gaiser [Ref. 10]. Component temperature measurements both without and with a 
shroud parallel to the heated surface were reported and experimental correlations were 
developed for different power levels and shroud spacings. A significant increase in tem- 
perature due to reduction in fluid velocity for the smallest channel spacing was observed. 

A numerical investigation of natural convection flow and heat transfer arising from 
a protruding heat source on a vertical plate was reported by Sathe and Joshi (Ref. 11]. 
Numerical predictions were obtained for a wide range of appropriate Ra and Pr numbers 
and substrate to fluid thermal conductivity ratios. Computed protrusion surface tem- 
peratures compared favorably with available experimental results from [Ref. 8]. A nu- 
merical investigation of two dimensional natural convection flow and heat transfer from 
a substrate-mounted flush heat source immersed in a liquid-filled square enclosure was 
conducted by Haukenes [Ref. 12]. Numerical predictions were in favorable agreement 


with the experimental data. 











While the transport characteristics during steady state have been investigated for 
both discrete flush and protruding heat sources, very little information is available dur- 
ing transient phases of operation. In many applications power cycling mav take place. 
The present investigation examunes both numerically as well as experimentally the tran- 
sient response of discrete heat sources flush mounted on a substrate. The experiments 
were conducted with the same test surface as of Gaiser [Ref. 10], and comparison be- 
tween numerical model and experiments were made for the heat source closest to the 
leading edge of the surface. 

The experimental studies reported here examine changes in the natural convection 
flow and heat transfer characteristics of a three column arrav of 45 flush mounted 
heaters on a vertical test surface. The power to the bottom heater of the center column 
oscillated with time while power to all the other heaters was kept constant. Component 
surface temperatures were measured over a period of time for a range of average power 
inputs from 0.2 W - 1.0 W. Various amplitudes and frequencies of the periodic variation 
were examined. Also reported here are the effects of the presence of a shroud. Numer- 
ical computations were performed to study heat transfer from a single heat source. 


Specific goals of this study were: 


e To measure oscillations in heater surface temperatures in response to input power 
fluctuations for a range of average component power levels, amplitudes and fre- 


quencies of power pulsations. 


e To examine the same effects in the presence of a parallel shroud placed next to the 


test surface. 


¢ To examine the effects of power pulsations in the bottom component on the 


steadily powered downstream components. 





¢ To numerically investigate the two-dimensional transient natural convection flow 
and heat transfer for a single flush heat source and compare the results with the 


experimental measurements of component temperatures. 


This study advances the work reported by Gaiser (Ref. 10] and Haukenes [Ref. 12] 


and is part of the study of the natural convection liquid immersion cooling of electronic 


equipment. 





II. EXPERIMENTAL SETUP AND PROCEDURE 


Detailed descriptions of the original experimental equipment design and con- 
struction can be found in Gaiser [Ref. 10] and further additions and changes done in the 
experimental setup in Haukenes [Ref. 12]. Only a brief review is presented and new 
noteworthy additions and changes are noted here. 

As seen in Figure 1, the assembly consists of a vertical test surface with a total of 
45 flush mounted rectangular heat sources in three columns of 15 elements per column. 
The test surface is immersed in a one cubic meter, plate glass tank filled with deionized 
water. A parallel shrouding plate, placed in front of the test surface to simulate a typical 
circuit board assembly, forms a channel with the test surface. Each heater has a 
copper-constantan thermocouple attached to its back. All the thermocouples are refer- 
enced individually to an ice bath as seen in Figure 2. Three thermocouples monitor the 
tank temperature at the top, middle and bottom. A Hewlett Packard Model 3852S data 
acquisition system measures the ice referenced voltage of each thermocouple to calculate 
temperatures and voltages delivered from the source and supplied to the heater to cal- 
culate power. The voltages are converted into temperatures and component input 
powers monitored by a Hewlett Packard Model 9153 computer interfaced to the data 
acquisition system. A filtration system consisting of four cartridges is used to ensure the 


purity and the resistivity of at least 0.7 megaohm.cm of the water inside the tank. 


A. MODIFICATIONS IN APPARATUS 
The most significant new feature was the design of a power assembly to provide 
pulsating input powers to the heating elements. Figure 2 shows a sketch of this assem- 


bly. All the heater strips in the left and right columns are powered by 








— 
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SECTION AA 


Heater 


Thermocouple Pad ae 


Figure 1. Test Assembly (After Gaiser [Ref. 10]). All Dimensions are in cm. 





a 0-40 V, 0-50 A HP 6269B D.C. power supply unit through four distribution panels. 
In the central column all heater strips except the one at the bottom are powered by a 
0-10 V, 0-10 A HP 6282A D.C. power supply unit through two distribution panels. 
There is a 2.01 ohm precision resistor in series with each heater resistor for current 
measurements. Details of the power distribuuon assembly are provided by Gaiser {Ref. 
10]. 

A HP-85 computer, HP-59501B power supply programmer and 0-40 V, 0-1.5 A 
HP-6289A D.C. power supply units were added to the experimental setup to generate a 
periodic power input. The transient power input is provided only to the heater at the 
bottom of the central column. Appropriate programs were developed for the HP-85 
computer to provide transient power to the bottom heater. 

The following relation was used to determine the power: 


(As seen from Figure 3) 








where the symbols are defined in the nomenclature. See Appendices B-F for program 
listings for power input and data reduction, developed for both steady and transient 


Cases. 
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Figure 3. Power Distribution to The Heater at The Bottom of The Centeral Column 











B. EXPERIMENTAL PROCEDURE 

The preperation for each set of experiments started by stirring the water in the 
immersion tank with a propeller tvpe stirrer for 5 to ‘0 minutes to achieve uniform 
temperature distribution. Crushed ice was placed into the thermocouple reference bath. 
A mercury thermometer was also inserted into the ice bath to make sure that the bath 
was at 0 °C +0.05 °C . The data acquisition system was turned on and warmed up. 
Typically in about one hour th: data acquision system warmed-up and also the 
immersion tank was quiescent. The temperature acquisition program (see Appendix D) 
was loaded and the temperatures of the heater surfaces and immersion tank were meas- 
ured, before the power was supplied to the heaters. Typically all heaters and the tank 
were within 0.05 °C of each other prior to power input. 

The experimental study consisted of three major parts. The heater surface temper- 
ature measurements for three different cases were made as presented in Table I. In the 
first case thermocouple measurements were taken every fifteen minutes until successive 
measurements remained unchanged. Once steady state was achieved, temperature 
measurements and heater voltage readings were recorded and saved for further process- 
ing. The computer program to generate a prescribed periodic power input, contained in 
Appendix B was then loaded and run. Once a steady periodic response was evident, 
temperature oscillations with time were sampled over a period of forty seconds with a 
sampling frequency of 180/40 Hz. The inpui po.cr pattern was also recorded over a 
period of seventy seconds with a sampling frequency of 350/70 Hz. 

As the amplitudes of the power variation increased, the accuracy in digitization of 
the input power by the data acquisition system decreased. While the power input to the 


heater was as programmed, the measurement of it by the HP-85 was sometimes slower 


than required. Another problem encountered was the monitoring of the voltage drops 








across both the power supply and the resistor at the same time in order for the calcu- 
lation of the power supplied to the heater. 

In the second catagory of experiments in Table 1, the power to all heaters was in- 
creased incrementally. The output from the power calculation program was checked, 
until the desired power level was reached. The power and temperature measurements 
were recorded when steady state was attained. The periodic input power program was 
then loaded and run. Besides monitoring the temperature oscillations with time on the 
bottom heater, downstream elements with uniform power were also sampled to detect 
any temperature oscillations with time, over a period of five minutes. 

In the third class of experiments in Table 1, the same procedure as for the second 
set was repeated. The temperatures of each heater element both for steady and transient 
cases were recorded. Additional runs for different average power levels and with differ- 
ent periods of applied power were also performed. For the first group of these runs, the 
effects of changing the pulsation frequency and for the second group, the effects of 
changing the amplitude of input power were examined. For the third group of runs, the 
effects of providing all the heaters with power different than supplid to the bottom one 


were studied. 
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Table 1. EXPERIMENTAL PROCEDURE FOR DIFFERENT CASES 


CASE PROCEDURE 


CASE 1}: Single 
heater powered 


CASE 2: Centeral 


column powered 


CASE 3: All heat- 


ers powered 


Constant power is ap- 
plied to the heater at 
the bottom of centeral 
column, until steady 
state is achieved. 


The power input is then 
modulated in time. 


Constant power is ap- 
plied to all heaters until 
steady state is achieved. 


The power to the heater 
at the bottom of cen- 
teral column is then 
timewise modulated 
while keeping the power 
to all other heaters 
steady. 


Same as CASE 2 
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POWER 
(Watts) 


049, .082, .2, 
.235, .4, .6, .8, 
1. 


Same as 


CASE 2 





SHROUD 
SPA.(mm) 


No shroud, 
3 


No shroud 


(steady and 


transient), 


3,6,9 (all 


steady) 


Table 2. EXPERIMENTAL PROCEDURE FOR SPECIAL RUNS 


CASE PROCEDURE POWER (Watts) 


Change in period of in- CASES I and 2 in Table 1 2, .6 















put power are repeated by changing the 





period from 20 sec. to 40 sec. 


and then 10 sec. 


Change in amplitude CASE | and 2 in Table | are 











of applied power repeated by changing the 





amplitude of applied power 


to selected values. 














Same as CASE 2 in Table 1 





Different power to .024 (bot)- .275 (others), 





bottom heater than the .049 (bot)- .45 (others), 





others in central col- .082 (bot)- 1.27 (others), 





umn .235 (bot)- 2.06 (others) 








If. EXPERIMENTAL RESULTS 


A. SINGLE HEATER ELEMENT DATA 

The presentation and discussion of the single heater element data during this study 
is divided into two major parts. The first part refers to the case where the heater element 
is powered with uniform power, until steady state condition is reached. Within this part, 
heater tempe.ature data from the thermocouple at the center of the heater surface are 
recorded for 40 sec. with a sampling frequency of 180/40 Hz. In the second set of ex- 
periments after an initial steady state condition is reached within the heater element, its 
input power is switched to a sawtooth form. In this form, input power is increased lin- 
early then decreased in the same fashion between maximum and mimumum amplitudes. 
Applied transient power has a period of 20 sec. The surface temperature of the heater 
element is recorded using the the same sampling frequency. 

These two sets of measurements are repeated for five average power levels, ranging 
from 0.2 W to 1.0 W. Measurements are made without shroud and also with a shroud 
at 3 mm spacing from the test surface, in order to determine the effects of power level 
and channel spacing on the heat transfer performance. In order to compare measured 
surface temperatures with the numerical computations, the heater center temperatures 
are presented as differences between the heater surface and the ambient. 

1. Measurements in The Absence of Shroud 

In Figures 4-8 the temperature variations with time are displayed for various 
average power inputs with a fixed period of 20 sec. In Figures 9-12 the effects of 
changing the period of applied power, for 0.2 W and 0.6 W average power levels are 
examined. In Figures 13-16 the effects of changing the amplitude of the applied power 


to 0.4 W and 0.8 W are presented. From Figures 4-12, the trend of increasing average 
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surface temperature with increasing average power is as expected. The period of surface 
temperature variation, noted almost as 23 sec. from these figures, stays approximately 
constant for all power levels. It is also seen that the oscillations in surface temperature 
increase with an increase in both amplitude and the period of the applied power vari- 
ation. For example in Figures 13 and 14 for an average input power of 0.4 W when the 
power variation is within 0.55 W and 0.29 W, the maximum difference in AT is 0.8 °C. 
This difference decreases to 0.25 °C, if the power variation is kept within 0.43 W and 
0.36 W. For 0.6 W average power level in Figures 15-16 when the period of applied 
power is increased to 40 sec., the difference in AJ increases to 1.0 °C. The difference in 
AT decreases to 0.3 °C when period is reduced to 10 sec. For higher power levels, the 
effect of amplitude and period of applied power is more evident. 
2. Measurements With Shroud 

Figures 17 and 18 show AT versus time at mean power levels of 0.2 W and 1.0 
W for shroud spacing of 3 mm. When compared to the unshrouded measurements at 
the same power levels [Figures 4 and 8], it is seen that the AT values for a single element 


remain almost unchanged. 
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Figure 4. Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations at 0.2 W. Pulsation Amplitude and Period are 0.06 W 


and 20 sec 
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Figure 5. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations at 0.4 W. Pulsation Amplitude and Period are 0.07 W 


and 20 sec 
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Figure 6. 
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Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations at 0.6 W. Pulsation Amplitude and Period are 0.1 W 


and 20 sec. 
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Figure 7. Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations at 0.8 W. Pulsation Amplitude and Period are 0.12 W 


and 20 sec 
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Figure 8. Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations at 1.0 W. Pulsation Amplitude and Period are 0.13 W 


and 20 sec 
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Figure 9. Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations with a Period of 10 sec. at 0.2 W. Pulsation Amplitude 


is 0.06 W. 
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Figure 10. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations with a Period of 40 sec. at 0.2 W. Pulsation Amplitude 


is 0.06W. 
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Figure 11. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations with a Period of 10 sec. at 0.6 W. Pulsation Amplitude 


is 0.1 W. 
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Figure 12. 
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Figure 13. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations with Mean Power of 0.4 W, Amplitude of 0.03 W and 


Period of 20 sec. 
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Figure 14. 
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Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations with Mean Power of 0.4 W, Amplitude of 0.13 W and 


Period of 20 sec. 
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Figure 15. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations with Mean Power of 0.8 W, Amplitude of 0.06 W and 
Period of 20 sec. 
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Figure 16. Oscillations in Surface Temperature of Buttom Heater Due To 'nput 
Power Pulsations with Mean Power of 0.8 W, Amplitude of 0.19 W and 
Period of 20 sec. 
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Figure 17. Oscillations in Surface Temperature of Bottom Heater Due To Input 
Power Pulsations with Mean Power of 0.2 W, Amplitude of 0.06 W and 


Period of 20 sec. 
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Figure 18. Oscillations in Surface Temperature of Bottom Heater Due To Input 


Power Pulsations with Mean Power of 1.0 W, Amplitude of 0.13 W and 
Period of 20 sec. 


29 








B. SINGLE COLUMN DATA 


In another set of experiments the effects of transient power, applied to the heater 
at the bottom of the central column were examined. The remaining fourteen heaters in 
the column were steadily powered. The main point of interest was to look for changes 
in the heat transfer characteristics compared with those existing with steady power to 
all heaters. The same experimental procedure that is followed in the previous section is 
used for powering the bottom heater. The steady and transient temperature data were 
taken for five power levels, ranging from 0.2 W to 1.0 W both with and without a 
shroud. The surface temperatures of heaters Ist, 10th and 13th were monitored once 
steady power of bottom heater was switched to the periodic form. Temperature meas- 
urements were made with a sampling frequency of 300/335 Hz. over a period of 335 
seconds. 

1. Measurements In The Absence of Shroud 

Figures 19-23 show AT versus x/L at steady input power levels of 0.2 W, 0.4 

W, 0.6 W, 0.8 W and 1.0 W respectively. From the plots of data, the trend of increasing 
surface temperature along the test surface is readily apparent. It is also seen that at a 
given location, surface temperature increases with rising power level. 

The steady temperature levels were used to determine the non-dimensional heat 
transfer characteristics. Rayleigh and Nusselt numbers were calculated as: 

Ra = (gf@P) / (akwv), Nu = hl/4,. Figures 24 and 25 are plots of Nu versus Ra 
for 0.2 W and 0.8 W respectively. The local values of h indicate a decreasing trend 
moving from the bottom heater to the top. It is also obvious that h values increase with 
rising power level at a given location due to the stronger flow. 

The transient temperature data for the bottom heater over a period of 40 sec. 


were the same as for the single heater element examined in the previous section. Figures 
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26 and 27 are the temperatures of Ist, [Oth and [3th heaters, over a period of 335 sec. 
for 0.4 W and 0.8 W input power levels respectively. As seen from the figures, oscillation 
in surface temperatures of downstream heaters, when power to the bottom heater is os- 
cillated with time are barely noticeable. 
2. Measurements With Shroud 

In the previous section, it was found that the effects of the shroud on the heat 
transfer from the single heater element were not apparent. In this section, the effect of 
shroud is investigated with the entire central column powered. Surface temperature 
measurements of heaters were recorded at five different power levels for 3 mm shroud 
spacing. The results are compared with the unshrouded measurements. 

The temperature rise above ambient for each heater at average power levels of 
0.2 W and 0.8 W for 3 mm shroud spacing are seen in Figures 28 and 29. From these 
plots temperature values higher than the unshrouded case measurements are seen be- 
yond the sixth heater element from the bottom. When the power levels increase, differ- 
ences in AT values also increase compared to the unshrouded case. These measurements 
are in agreement with Knight [Ref. 13] and Gaiser (Ref. 10]. While the thermal layer 
near the bottom component is thin enough not to be effected by the presence of the 


shroud, the same is not true further downstream. 
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26 and 27 are the temperatures of Ist, 10th and 13th heaters, over a period of 335 sec. 
for 0.4 W and 0.8 W input power levels respectively. As seen from the figures, oscillation 
in surface temperatures of downstream heaters, when power to the bottom heater is os- 
cillated with time are barely noticeable. 
2. Measurements With Shroud 

In the previous section, it was found that the effects of the shroud on the heat 
transfer from the single heater element were not apparent. In this section, the effect of 
shroud is investigated with the entire central column powered. Surface temperature 
measurements of heaters were recorded at five different power levels for 3 mm shroud 
spacing. The results are compared with the unshrouded measurements. 

The temperature rise above ambient for each heater at average power levels of 
0.2 W and 0.8 W for 3 mm shroud spacing are seen in Figures 28 and 29. From these 
plots temperature values higher than the unshrouded case measurements are seen be- 
yond the sixth heater element from the bottom. When the power levels increase, differ- 
ences in AT values also increase compared to the unshrouded case. These measurements 
are in agreement with Knight [Ref. 13] and Gaiser (Ref. 10]. While the thermal laver 
near the bottom component is thin enough not to be effected by the presence of the 


shroud, the same is not true further downstream. 
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Figure 21. x/L versus AT at 0.6 W Input Power 
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Figure 22. x/L versus AT at 0.8 W Input Power 
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Figure 26. 


Figure 27. 
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Monitoring The Surface Temperatures of Selected Heaters at 0.8 W 
When Input Power To The Bottom Heater Oscillated. Oscillation Am- 


plitude and Period are 0.07 W and 20 sec. respectively. 
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Figure 28. x/L versus AT at 0.4 W With 3 mm Shroud Spacing 
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Figure 29. x/L versus AT at 0.8 W With 3 mm Shroud Spacing 
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C. DATA WITH ALL THREE COLUMNS POWERED 

A third set of experiments was undertaken to examine the effects of placing flanking 
columns of heaters on each side of the central column. The three column results are 
discussed next. In addition to the measurements for the ceniral coluni, the surface 
temperatures of top and bottom heaters of both side columns were also monitored for 
335 sec. with a sampling frequency of 300/335 Hz. These transient measurements were 
repeated for average power levels ranging from 0.2 W to 1.0 W with 3mm shroud spacing 
and also without shroud. Additional runs were performed with steady power inputs, 
within the same power range and shroud spacings of 6 and 9 mm. 

1. Measurements In The Absence of Shroud 

The surface temperatures in each of the three columns are plotted to indicate 
changes in convective heat transfer characteristics as compared with the single column 
case. Figures ’9 and 31 are AT versus x/L plots at steady power levels of 0.2 W and 0.8 
W respectively. from these plots it is seen that at any x/L the central column has 
slightly higher AT values when compared to the right and left columns. Referring back 
to Figures 19 and 23 for a single column, it is apparent that the AT valies for central 
column components with all three columns powered are similar to the values obtained 
for single column powered case. Comparing the heat transfer coefficients for chese ex- 
periments with the single column values it was seen that the three column data produced 
almost the same values of h. 

As was done with the single column, the surface temperature of the bottom 
heater of the central column was monitored for 40 sec. with a sampling frequency of 
180/40 Hz, after modulating its input power with a time period of 20 sec. Once again 
the temperatures of the three heaters above it and the heater at the top of each column 
were also monitored for 335 sec. with a sampling frequency of 300/335 Hz. Figures 


32-34 are plots of AT versus time for mean power 0.2 W, amplitude 0.06 W and Figures 
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35-37 for average power 0.8 W, amplitude 0.12 W. From the figures, no oscillation in 
the surface temperature in any heater indicates the damping of the effect of the input 
power variation in the bottom heater is observed. It is also noted that the surface tem- 
perature oscillation of bottom heater due to periodic power, follows the same pattern 
as for a single heater. 
2. Measurements With Shroud 

Measured heater surface temperatures are presented for the central column in 
Figures 38-40 at various shroud spacings and power levels. It is seen that the shroud 
has a significant effect on the heater temperatures only for small spacings. For spacings 
greater than 6 mm, the heater temperatures are almost equal to the unshrouded test 
case. Figure 41 compares the surface temperatures for 3 mm shroud spacing at various 
power levels. Figure 42 compares the measured surface temperatures without shroud 
with the values for 3 mm shroud spacing at 0.4 W and 0.8 W power levels. From the 
figure it is seen that as the shroud spacing decreases, surface temperatures of central 
column heaters are increased. Figure 43 shows the surface temperatures of all heaters 
for 3 mm shroud spacing at 0.4 W. From the figure it is seen that for this shroud spac- 
ing, surface temperatures of the central column heaters are higher than the measured 
values for the side columns at the same x/L. Fluid entrainment with shroud in place is 
impeded in the direction normal to the test surface. Lateral fluid movement brings 
warmer fluid from the flanking columns towards the central column. This results in an 


increase in temperature of the central column heaters. 
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Figure 30. x/L versus AT at 0.2 W Input Power 
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Figure 31. x/L versus AT at 0.8 W Input Power 
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Figure 32. Monitoring The Surface Temperatures of Selected Heaters at 0.2 W 
When Input Power To The Bottom Heater Oscillated. Oscillation Am- 

plitude and Period are 0.06 W and 20 sec. respectively. 
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Figure 33. Monitoring The Surface Temperatures of Selected Heaters at 0.2 W 


When [Input Power To The Bottom Heater Oscillated. Oscillation Am- 


plitude and Period are 0.06 W and 20 sec. respectively. 
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Monitoring The Surface Temperatures of Selected Heaters at 0.2 W 


When Input Power To The Bottom Heater Oscillated. Oscillation Am- 
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Figure 35. Monitoring The Surface Temperatures of Selected Heaters at 0.8 W 
When Input Power To The Bottom Heater Oscillated. Oscillation Am- 


plitude and Period are 0.12 W and 20 sec. respectively. 
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Figure 36. Monitoring The Surface Temperatures of Selected Heaters at 0.8 W 


When Input Power To The Bottom Heater Oscillated. Oscillation Am- 


plitude and Period are 0.12 W and 20 sec. respectively. 
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Figure 37. 
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Monitoring The Surface Temperatures of Selected Heaters at 0.8 W 
When Input Power To The Bottom Heater Oscillated. Oscillation Am- 


plitude and Period are 0.12 W and 20 sec. respectively. 
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x/L versus AT for Different Shroud Spacings At 0.2 W 
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Figure 41. x/L versus AT for 3 mm Shroud Spacing At Various Power Levels 
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Figure 42. Comparison of Surface Temperatures Measured Without And With A 
Shroud at 3 mm For 0.4 W and 0.8 W 
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Figure 43. Effect of Shroud on Surface Temperatures of Heaters on Different Col- 


umns of Test Surface with The Same x/L 
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IV. NUMERICAL MODEL AND COMPARISON WITH 
MEASUREMENTS 


Based on the experimental results presented in the previous chapter it is clear that 
in the absence of shroud the lateral fluid entrainment is not very strong. With only a 
single component powered, the transport away from the ends can then be modelled as 
two dimensional. In the following, a computational model of natural convection trans- 


port from a single flush heater is presented. 


A. MATHEMATICAL FORMULATION 

The schematic sketch of the configuration numerically investigated is shown in Fig- 
ure 44. The flow is assumed to be two dimensional, resulting from a heat source, flush 
mounted on a vertical substrate and located at the middle of the surface. The shaded 
region represents the substrate, and the darker region in it represents the heater. The 
heater, the substrate and the fluid have constant thermophysical properties. The vertical 
coordinate x is measured upward from the base of substrate and u is the velocity com- 
ponent in this direction. The y coordinate is measured horizontally from the surface and 
Vv is the velocity component in the same direction. Isothermal impermeable enclosure 
walls with no-slip i.e u= v= T=0 are prescribed as boundary conditions. The governing 
equations assuming no viscous dissipation and pressure energy terms and the Boussinesq 


approximations to be true, are: 


FLUID REGION : 
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Figure 44. Numerical Model (After Haukeness [Ref. 10]) 
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B. COMPUTATIONAL PROCEDURE 

The information contained in the differential equations is discretized which makes 
it possible to replace the governing partial differential equations with algebraic 
equations, whick can be solved with relative ease. A finite difference scheme is used for 
discretizing the equations. The calculation domain is divided into a number of control 
volumes such that there is one control volume surrounding each grid point. The tem- 
perature and pressure components are located at the main grid points. The location of 
the velocity components is staggered and defined on the control volume faces. The 
power law profile that provides a very good fit to the exact expression, is used for the 
spatial variation of the dependent variables between grid points. The harmonic mean 
formulation that provides physically meaningful results even with large changes in ma- 
terial properties is followed for the calculation of the interface diffusivities. The details 
of the discretization process can be found in Patankar [Ref. 14]. 

The discretized, nominally linear equations are solved by line-by-line TDMA ( Tri- 
Diagonal Matrix Algorithm ) method in which the coefficients are iteratively recalcu- 
lated and the process repeated until convergence. ° Convergence is based on a balance 
of the rate of energy generated in the heater and the rate of energy leaving the enclosure 
walls. The SIMPLER (SIMPLE-Revised) algorithm that employs two different 
equations: one for the task of correcting the velocities, the other for obtaining a new 
pressure field, is followed during solution as outlined in [Ref. 14]. Even though separate 
equations are written for solid and fluid regions, numerically the solid region in the mo- 
mentum equation is simulated by specifying a very large viscosity. The heat source term 
and property values in the energy equation are specified appropriately in the various 
regions. As a result, only one momentum and one energy equation are so!ved through- 


out the computational domain. 
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Numerical computations were performed on 48x56 control volumes. Different time 
steps were tried ranging from 0.) sec. to 5.0 sec. to determine their effect on the scl- 
utions. Calculations for the modulated power input were started with u, v and T values, 
determined previously for a steady case with the same average power. The initial profiles 
were obtained by running the transient code with a very large time step. 

The compuiations were performed initially on a IBM-370/3033 computer and later 
on an Amdahl Computer. Appendix | shows the problem dependent portion of the fi- 


nite difference program used to obtain the the numerical results. 
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C. COMPARISON OF COMPUTATIONS WITH EXPERIMENTAL 
MEASUREMENTS 

Numerical investigation of transient heat transfer consisted of two parts. For the 
first, initial transient response due to sudden step power input was studied. Figures 45 
and 46 show representative transient flow and temperature contour developments for 
8.69 W'm step input. We note that for the heater size used in the experiments this 
power input corresponds to 0.2 W. The figures from left to right show the contours at 
10, 44, 80 and 120 seconds following start of heating. At 10 sec. the isotherms are al- 
most parellel to the heater indicating a diffusion dominated response. A single clockwise 
cell exists next to the heater. In the second and third figures, the secondary cell appears 
and grows. Flow becomes relatively faster and spreads further away. Also at 120 sec. 
the flow begins to resemble the fully developed condition seen in Figure 54 for the same 
input power. 

For the second part of numerical study, the transient response following a periodic 
input power into the heater was studied. The amplitude and time period of imposed 
power oscillations were varied. Time steps from 0.1 sec. to 5.0 sec. were investigated. 
Surface temperatures calculated for different time steps are presented in Table 3. Tran- 
sient variations of surface temperatures for different time steps are displayed in Figures 
47 - 50. The surface temperature patterns are similar for all tume steps. The highest 
difference in calculated temperature values was 2.25 %, when 5.0 sec. time step was used 
instead of 0.1 sec. Difference in calculated surface temperatures was smaller than 0.45 
%o, if calculations were performed with 0.5 sec. time step rather than 0.1 sec. In order 
to optimize the computation times without compromising significantly on the accuracv, 
all numerical calculations were carried out with a time step of 0.5 sec. 

Numerical computations were also performed for two different sets of heating ele- 


ment properties since the actual element properties in the experiments were not easily 
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quantifiable. In the experiment, the composite heater element consisted of an inconel 
ribbon attached to a kapton substrate. In order to cover the range of property vari- 
ations between the two materials, the computations were repeated, with the properties 
of incone] and of plastic separately for the heater [Ref. 15]. Figures 51 and 52 show the 
numerically calculated surface temperatures with time for average power levels of 0.2 
W and 1.0 W for plastic heater and Figure 53 for inconel heater at 1.0 W average power 
level. While the mean temperatures in Figure 53 are lower than in Figure 52, the am- 


plitudes of temperature oscillations remain approximately the same. 


Table 3. EFFECT OF SELECTING DIFFERENT TIME STEPS ON COM- 
PUTED QUANTITIES 
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Table 4. THERMOPHYSICAL PROPERTIES OF MATERIALS IN HEATER 
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Figure 45. Development of Flow During Initial Transient Due To 0.2 W Step Input 
Power. The Figures are at (i)10, (ii)44, (iii)80 and (iv)120 seconds. 
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Figure 46. Development of Temperature Contours Due To 0.2 W Step Input 


Power. The Figures are at (i) 10, (ii)44, (iii)80 and (iv)120 seconds. 
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Surface Temperature Calculated with Time Step of 0.1 sec. 


Figure 47. 
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Figure 48. Surface Temperature Calculated with Time Step of 0.5 sec. 
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Figure 49. Surface Temperature Calculated with Time Step of 1.0 sec. 
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Figure 50. 


Figure 51. 
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Surface Temperature Calculated with Time Step of 5.0 sec. 
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Figure 52. 


Figure 53. 
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COMPARISON 

The bottom heater in the central column (Section AA in Fig.4) was chosen for a 
comparison with numerical results. For 0.2 W average power input with 0.06 W ampli- 
tude and 20 sec. pulsations, the steadv periodic temverature response, calculated at the 
back of the model heater, was 27 % lower than the measured for an inconel heater and 
38 % higher than the measured value for a plastic heater. The oscillation amplitudes in 
calculated surface temperatures agreed within approximately 33 % with the exper- 
imentally measured values. For the 1.0 W average power case, the calculated surface 
temperature of an inconel heater was 13 % lower than the experimentally measured 
value. The calculated value for a plastic heater was 51 % higher than the measured one. 
The amplitudes of oscillations in calculated surface temperatures agreed within 33 % 
with the experimental values. 

The percent differences, presented in the following between numerically calculated 
and experimentally measured surface temperatures, arc for the case when the heater 
material is assumed to be inconel. Figures 54 and 55 show the corresponding temper- 
ature contours and streamlines for the average power input levels of 0.2 W and 1. W. 
The effect of change in the period of the input power pattern was also examined nu- 
merically. Figures 56 and 57 show the surface temperature variation with time for 0.2 
W average power with different periods. The trend for increasing period of the measured 
temperature response due to increase in period of applied power is evident. When 
compared to the experimentally measured values, the calculated temperature oscillation 
amplitudes were almost 33 % lower for all power levels. Figures 58-60 show the effects 
of the period of applied power for a higher average power level of 0.6 W. 

The effect of change in the amplitude of the applied power variation was also in- 


vestigated. Computations revealed increasing differences between peak values of calcu- 
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lated heater surface temperatures due to increase in amplitude of applied power. The 
calculated surface temperatures versus time for 0.4 W and 0.8 W average power levels 
with changing amplitudes are presented in Figures 61 and 62. The corresponding tem- 
perature contours and streamlines are plotted in Figures 63 and 64. The calculated 
steady mean temperature was almost 21 °%o lower for 0.4 W and 13 % lower for 0.8 W 
compared to the experimentally measured values. The calculated oscillation amplitudes 
in surface temperatures agreed only within 67 % with the experimental values for both 
power levels. 

This difference can be ascribed to the larger thickness of the heater model than the 
actual value. The actual thickness was about .195 mm. However a very fine grid will 
be required if this size is to be taken in the model. As the thickness of the heater in the 
computational model was reduced, the oscillation amplitudes increase. This trend is 
presented in Figure 66 with the depth of the heater reduced to half of the previous value. 
If the actual properties of heater element which are combination of inconel and kapton 
properties were available, it is obvious that the agreement with the experiment will im- 
prove. The experimentally measured temperature values are bracketed by the compu- 
tations, as expected. 

Additional computations were performed with k values between inconel and plastic 
to achieve better agreement between the model and experiments. Figure 65 shows the 
surface temperatures calculated with k, = 1.45 W/m.°C. It is seen that calculated val- 
ues are reasonably close to the experimental values. A particular run was performed 
with k= 1.0 W/m-K and heater thickness of .195 mm. Computed results are presented 
in Figure 67. For these choices of thermophysical properties and component dimensions 


close agreement is seen with the experimental data in Figure 4. 
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Figure 54. Temperature Contours and Streamlines with An Inconel Heater for 


Mean Power of 0.2 W, Amplitude of 0.06 W, Period of 20 sec at 15th sec 
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Temperature Contours and Streamlines with An Inconel Heater for 
Mean Power of 1.0 W, Amplitude of 0.13 W, Period of 20 sec at 15th sec 


Figure 55. 
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Figure 56. 


for 0.2 W Average Power and Amplitude of 0.06 W. 
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Surface Temperature Due to Power Pulsation with a Period of 10 sec. 


Figure 57. 


for 0.2 W Average Power and Amplitude of 0.06 W 
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Figure 58. Surface Temperature Due to Power Pulsation with a Period of 40 sec. 


for 0.6 W Average Power and Amplitude of 0.1 W 
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Figure $9. Surface Temperature Due to Power Pulsation with a Period of 20 sec. 


for 0.6 W Average Power and Amplitude of 0.1 W 
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Figure 60. Surface Temperature Due to Power Pulsation with a Period of 10 sec. 


for 0.6 W Average Power and Amplitude of 0.1 W 
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Surface Temperature Due to Power Pulsation with Different Amplitudes 
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Figure 62. 


for 0.8 W Average Power 
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Temperature Contours and Streamlines for 0.4 W Average Power, 0.07 


Figure 63. 


W Amplitude and 20 sec. Period of Oscillation 
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Figure 64. Temperature Contours and Streamlines for 0.8 W Average Power, 0.12 


W Amplitude and 20 sec. Period of Oscillation 
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Figure 65. Effects of Thermophysical Properties on Calculated Temperatures 
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Figure 66. Effects of Heater Dimensions on Calculated Temperatures 
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Figure 67. 
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V. CONCLUSIONS 


Following conclusions are drawn based on this experimental and numerical study 
about the heat transfer effects of input power pulsations. 

The experiments: 

The oscillation amplitudes in surface temperature increase with an increase in both 
the amplitude and the time period of the applied power. The period of temperature os- 
cillation remains almost the same for all shroud spacings at a particular power level. It 
was observed that the effects of periodic power, applied to the heater at the bottom of 
the central column are barely noticeable on the heat transfer characteristics of the other 
steadily powered heaters in the column. There is no noticeable oscillation in surface 
temperatures of downstream heaters when power to the bottom heater oscillates in time. 
With all three columns powered under steady conditions and with a shroud in place, the 
central column heaters at any x/L have slightly higher AT values when compared to the 
right and left columns. For spacings greater than 6 mm, the heater temperatures are 
almost equal to the unshrouded test case. These conclusions regarding the steady re- 
sponse are in agreement with Knight [Ref. 13] and Gaiser [Ref. 10]. 

The computations: 

Numerical computations of the transient response due to both a sudden step power 
input and a periodic power input were conducted. Transient flow and temperature pro- 
files were obtained. Computations were performed for two different heating element 
properties since the actual element in the experiments was a composite and its properties 
were not easily quantifiable. Thermophysical properties of the heating element and its 
thickness had a pronounced effect on calculated temperatures. As the thickness of the 


heater element was reduced, the temperature oscillation amplitudes increased. Once the 
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variations in the thermophysical properties and appropriate thickness of the heating el- 


ement Were taken into account in the numerical study, ihe predicted temperatures were 


in reasonable agreement with measurements. 








VI. RECOMMENDATIONS 


In continuation of this study it is recommended that the following aspects be further 


investigated: 


¢ The effects of pulsating power inputs to more than one heater on the heat transfer 


from the other steadily powered heaters should be considered. 


e Investigate heat transfer enhancement due to possible turbulence effects caused bv 


high amplitude power pulsations at relatively high average power levels. 
e Use a multiple heater numerical model to determine heater interaction effects. 


e Analyze transient conditions for both one and multiple heater numerical models in 


the presence of a shroud. 


¢ Numerically study the effects of thermophysical properties of substrate on the am- 


plitudes of oscillations of surface temperature. 
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APPENDIX A. SUBROUTINE USER 


SUBROUTINE USER 
LOGICAL LSOLVE,LPRINT,LBLK,LSTOP,LSWICH 
COMMON F(66,66,10),P(66,66),RHO(66,66),GAM(66,66),CON(66,66), 
1 AIP(66,66),AIMC66,66),AJP(66,66),AJM(66,66),AP(66,66), 
2 X066),XU066),XDIF(66),XCV(66),XCVS(66), 
3 Y¥(66),YVC66),YDIF(66),YCV(66),YCVS(66), 
G YCVR(66),YCVRS(66),ARX(66),ARXJ(66),ARXJP(66), 
5 R¢666),RMNC66),SX(66),SXMNC66),XCVIC66),XCVIP(66) 
COMMON DU(66,66),DV(66,66), FV(66),FVP(66), 
1 FX(66),FXMC66),FY(66),FYM(66),PT(66),QT(66) 
COMMON/INDX/NF,NFMAX,NP,NRHO,NGAM,L1,L2,L3,M1,M2,M3, 
LIST, JST, ITER, LAST, RELAX(13),TIME,DT,XL,YL, 
ZIPREF, JPREF,LSOLVE(10),LPRINTC13),LBLK(10),MODE,NTIMES(10),RHOCON 
COMMON/CHTL/LSTOP,LSNICH 
COMMON/SORC/SMAX,SSUM 
COMMON/COEF/FLOW, DIFF, ACOF 
COMMON’ TRANS/CTIME, DELT, TLAST 
COMMON/OLD/ VOLD(66,66),VOLD(66,66),E0LD(66,66), IBODY(66,66) 
DIMENSION UC66,66),V(66,66),PC(66,66),COOL(66,66) 
DIMENSION 11(66,66),EN(66,66) 
EQUIVALENCECFC(1,1,1),UC1,1)),CFC1,1,2),V01,1)),CFC1,1, a7 
DIMENSION TH(66),THUC66),THDIF(66),THCV(66),THCVS(66), 
EQUIVALENCE(X, TH), (XU, THU), CXDIF, THDIF), (XCV, THCV) 
1 = ,CXCVS,THCVS),CXL,THL) 
(C38 30 36 3 DE 3 36 3 3 3 3 36 36 IE IE 36 3E 38 3 30 3G 3K 3 IE IE BE 36 30 3 DE DE 38 3 3E BE DE 3 36 36 30 9 9 DE OE IE 3E 3E 3 3 3 3 3 3 IE DE DE 3 3 3 DE IE 3E 3E 36 3E 3 3 3K 
DIMENSION E(66,66),7(66,66),COND(66,66), 
$TOLD(66,66),RHOOLD( 66,66) 
EQUIVALENCE (FC1,1,4),E(1,1)) 
EQUIVALENCE (F(1,1,5),T(1,1)) 
EQUIVALENCE (FC(1,1,10),COND(1,1)) 


ENTRY GRID 

ULLLLL“ZZLZ = GRID PARAMETERS % 
=0.0078 

DTH=0.023 


xX 
x 
XN 
x 
xn 
x 
xX 
x 
N 
x 
x 
x 


AL1=0.5*H1 
AL2=3.%H1 
AL3=1.xH1 
ALG=0.5xH1 
AL5=4.5*H1 
AL6=0.5x*H1 
BL1=2.xH1 
BL2=2.*H] 
BL3=1.*H1 
BLG=2.*H1 
BL5=3.*H1 
DX1=ALI1/1N1X 
DX2=AL2/N2X 
DX3=AL3/N3X 
DX4=AL4/NG4X 
DX5=AL5/N5X 
DX6=AL6/N6X 
DY1=BLI/N1Y 
DY2=BL27H2Y 
DY3=BL3/N3Y 
DYG=BLG/NGY 
DYS=BLS/N5Y 
xUC2)=0.0 
LFIRST=3 
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10 


20 


30 


40 


50 


60 


11 


21 


31 


41 


51 


LLAST=N1X+2 

DO 10 I=LFIRST,LLAST 
XUCI) =XUCI-1)+DX1 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST+N2X 

DO 20 I=LFIRST,LLAST 
XUCT) =XUCI-1)4DX2 
CONTINUE 
LFIRST=LLAST+1 
LLASTSLLAST+NSX 

DO 30 I=LFIRST,LLAST 
XUCI) =XUCI-1)4+DX3 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST+N4X 

DO 40 I=LFIRST,LLAST 
XUCT ) =XUCI-1)+DX4 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST#N5X 

DO 50 I=LFIRST,LLAST 
XUCI) =XUCI-1)+DX5 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST+N6X 

pO 60 I=LFIRST,LLAST 
XUCI) =XUCI-1)+DX6 
CONTINUE 

LL=LLAST 
XL=XUCLLAST) 
yv(2)=0. 

LFIRST=3 

LLAST=N1Y+2 

DO 11 I=LFIRST,LLAST 
YVCI)SYVCI-1)4DY1 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST#N2Y 

DO 21 I=LFIRST,LLAST 
YVCI)=YVCI-1)4+D¥2 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST+NSY 

DO 31 I=LFIRST,LLAST 
YVCI)=YVCI-1)4DY3 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST+NGY 

DO G1 I=LFIRST,LLAST 
YVCI) FYVCI-1)+DY4 
CONTINUE 
LFIRST=LLAST+1 
LLAST=LLAST#NSY 

DO 51 I=LFIRST,LLAST 
YVCI)2YVCI-1)4+DY5 
CONTINUE 

M1=LLAST 
YLEYVCLLAST) 

RETURN 


ENTRY START 
NHTRX=NIX4+N2X4N3X4+2 
NHTRY=N1LY4N2ZY+2 
NHTRX1=NHTRX+NSX-1 
NHTRYL=NHTRY+NSY~-1 
NBOSTX=NIX4N2X+2 
NBOFNX=N1X4+N2X+N3X4+NGX41 
NBOSTY2N1Y+2 
NBOFNY=N1Y+N2Y+N3SY4NGY41 


PRINT %,*INPUT LAST TIME VALUE IN SEC. AND* 


PRINT *,*TIME STEP INTERVAL' 
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READ ®, TLAST, DELT 
PERI=20.0 
DT=DELT 
CTIME=0.0 
DTIME=0.0 
PRINT ¥%, ‘INPUT DESIRED TIME FOR CONTOUR PLOTING IN SEC. ! 
READ *,TRINT 
REWIND (8) 
READ(8) TM 
D0 3000 I=1,L1 

3000 READ (C2) (UC1,J),J=1,M1) 
DO 3001 I=1,L1 

300i READ (8) a ue J=1,M1) 
DO 3002 I=1 

3002 READ (3) CEL, ay J=1 1) 
DO 3003 I= 

3003 READ (8) (FUL ed) 5), de st) 
DO 3004 I=1,L11 

3004 READ 2) (PCI,J),J=1,M1) 
DO 3005 I=1,L1 

3005 READ (8) (STOCI,J),J=1,M1) 
CP=4178., 
CPS=1500. 
CPHTR=439. 
RHOF=995. 
RHOS=1180. 
RHOH=8510. 

er9-807 


ALFAM=0.4 
DYNVIS=774.2E-6 
CONFLD=0.6027 
CONPLX=0.20 
CONHTR=11.7 
PT=3.14159256 
PR=DYNVISXCP/CONFLD 
VIRITECH,%) "PR=",PR 
RETURN 


aang 


ENTRY SETTIM 

PRINT %, "CTIME=', *DTIME=", DI IME 
IFCDTIME.LT.10.1) THEN 
EIN=(.364+0.07*¥DTIME/10.)/WIDTH 


ELSE 

EIN=(.430-CDTIME-10.)%.07710. “WIDTH 

ENDIF 

ENGEN=EIN/( DEPTHH 

POW=ENGENXWI DTHRDEP THRHL 

IPRTIM=0 

WRITE (1,%) CTIME, ECNHTRX+] , NHTRY+6),POW 

IFCCDTIME.Gi.9.95). AND. (DTIME.LT.10.05)) IPRTIM=1 

IFCDTIME.GT.CPERI-.1)) THEN 

DTIME=0. 

IPRTIM=1 

ENDIF 

IF CIPRTIM.EQ.1) THEN 

DSUM=0.0 

POVIER=0.0 

POHWER1=0.0 

DTOT=0.0 

DO 504 I=2,Lt2 
POWER=CONFLD¥CXUCT4+1)-XUCT) RECT, 29700 YVC03)-YV(02) 972.) 
POWER] =CONFL DRC XUCT4+1)~XUCI) RECI,M297CCYVOM1I-YVCM2) 72.) 
DSUM=POWER+POWER1+DSUM 

$04 CONTINUE 

DSUM1=0.0 

POWER2=0 

POWER3=0.0 
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DO 506 J=2,M2 
PONER2 = CONFLD¥CYV(J41)-~YVC J) 8EC2, S970 CXUC 3) -XUC2) 72.) 
POWERS=CONFLDECYV(J+1-YVCJ) )HECL2, 3970 (XUCLI)-XUCL2) 72.) 
DSUM1=DSUM1+POWERZ+POHERS 

§06 CONTINUE 

HRITEC#,%) ENERGY IN=', EIN 

WRITEC#,%) "ENERGY OUT=*,DSUM+DSUMI 

WRITEC%,%) "CURRENT TIM=',CTIME 

ENDIF 

FARK=TRINT-CTIME 

IF CTRINT.EQ.CTIME) THEN 

pO 999 I[=1,4 

LPRINT(I)=. TRUE. 

CALL PRINT 

LPRINTC(I)=. FALSE. 

999 CONTINUE 

NTIME=1 

WRITEC6,%) "START WRITING ON UNIT 9° 

WRITE (%,%) "ONCE YAZDIM' 

REWIND (3) 

WRITE (3) NTIME,M1,L1 


WRITE (3) CYVCJ),J=2,M1) 

WRITE (3) (xXUCI),I=2,11) 

WRITE (3) CYC), J=1,M1) 

WRITE (3) (XCI), 121,41) 

WRITE (3) CTIME 

DO 300 J=1,Ml 

300 WRITE (3) (VCI,J),T=1,L12) 

1 J=1,M1 


DO 30 

301 WRITE (3) CUCI,J),I1=1,11) 
DO 302 J=1,M1 

302 WRITE (3) CECI, J), 121,01) 


DO 303 J=1,M1 

303 WRITE 2 Par, J),I=1,L1) 
DO 305 J=1,™ 

305 WRITE ae iStoct, J),Ts1,L1) 
DO 304 J 

304 WRITE (3) TEL, J,3),T=1,L1) 


DO 306 J=1,M1 

306 WRITE raf CIBODY(I,J),3=1,L1) 
WRITEC(6,%) *CLOSE UNIT 3° 
CLOSE (3) 


881 CONTINUE 
LSTOP=.FALSE. 
LSWICH=.FALSE. 
ITER=0.0 
DIIME=DTIME+DELT 
RETURN 


ENTRY DENSE 
RETURN 


ENTRY BOUND 
RETURN 


ENTRY OUTPUT 
IFCITER.NE.O) GOTO 500 
VRITEC*, 501) 
501 FORMATC' ITER', 6X, "SMAX*,8X,* SSUM',7X, "VELCFRONT)', 
16X." TCHEATER) ' 6X, "CURRENT TIME') 
§00 CONTINUE 
WRITEC®,502) ITER, SMAX,SSUM, VCNBOFNX41,NHTRY+6 ) 
1, ECNHTRX4+1,NHTRY+6),CTIME 
502 FORMAT (16, 1P5E12.3) 
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RETURN 


ENTRY GAMSOR 
IF (NF.EQ.3) RETURN 
IF (NF.EQ.1.0R.NF.EQ.2) THEN 
DO 400 J=1,Ml 
po 400 I=1,L1 
RHO( IT, J) =RHOF 
GAM(I,J)=DYNVIS 
concr,J)=0.0 
IBODYCI,J)=1 
400 CONTINUE 
DO 401 I=NBOSTX, NBOFNX 
DO 401 J=NBOSTY,NBOFNY 
GAM(I,J)=1.E+20 
401 @ONTINUE 


ENDIF 

IFCNF.EQ.2) THEN 

DO 402 I1=2,L2 

DO 402 J=3,M2 
TEM=CFY( JECT, J)) + CFYMCJ)RECI,J-1)) 
CON) =RHOFXGHBETA®( TEM-TSTART) 
CONC I,J) =ALFAXCONI+ALFAMXSTOCI,J) 
STO(1I,J)=CON] 

402 CONTINUE 


ENDIF 
IF(NF.EQ.4) THEN 
DO 403 J=1,M1 
DO 403 J=1,L1 
RHO(T, J) =RHOFXCP 
GAM(I,J)=CONFLD 
CON(I,J)=0.0 
403 CONTINUE 
DO 404 I=NBOSTX, NBOFNX 
DO 404 J=NBOSTY,NBOFNY 
RHO(T,J)=RHOSXCPS 
GAMCI,J)=CONPLX 
404 CONTINUE 
DO 405 J=NHTRX,NHTRX) 


DO 405 J=NHTRY,NHTRY] 
RHOCI,J)=RHOHXCPHTR 
GAM(1I,J)=CONHTR 
CON(I,J)=ENGEN 

405 CONTINUE 
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APPENDIX B. STEADY STATE AND TRANSIENT POWER PROGRAMS 
FOR HP_85 


POWER PROGRAM FOR \ 
HP-S3SG81B PROGFAMMER | 


‘ 
1 
4 
t 
1 
4G) TO SUPPLY STEADY VOLT! 
' 
' 70 THE HEATER AT THE | 
1 
t 


eaTTOM OF CENT. CML. |! 


29 OLS "INPUT MAX. AND AVE. VO 
LTAGE” 

6 TIHPUT Wil 

@ R1L=V1-111%1499+2500 

QB QUTPLT 736 USING “#,40" ; Ri 
@ ENO 





5@ OIM E6193), V1¢199) 

66 DISP "INPUT MAX. VOL® 

76 DISP "AND MIN. VOL" 

80 INPUT R1i,R2 

96 DISF "INPUT PERIOD" 

100 INPUT P 

110 D2. @s. 

126 22P/C 

138 NSINT(Z> 

148 K2N7z 

158 V1ICK) ZINTCR1/254196)%18+2000 
168 Y1¢(8)=INT(R2/2541096)%10+2008 
170 ECK>2h1 

180 EC @)>=Rz 

198 L2k-1- 
266 FOR I=1 TOL 
210 ECT>e(R1-R2)>/-K+EC I-19 
220 ECN-1ZECT) 

23 Boe ENTSECT>7253100)210028 
240 VICN-Teavicl: 
250 NEXT I 
260 M2N-1 
278@ BEEP 

230 FOR I=0 TO M 

298 OUTPLT 706 USING "#,.40" ; VI 

CID 

306 WAIT 56 

310 NEXT I 

328 GOTO 276 

338 ENO 











APPENDIX C. TEMPERATURE AND POWER ACQUISITION ROGRAM 


FOR THE HEATER AT THE BOTTOM OF THE CENTERAL COLUMN 


12 
20 
38 
40 
50 
68 
70 
80 
90 
100 
110 
120 
130 
142 
15 
162 
178 
180 
19@ 
200 
210 
220 
230 
240 
250 
260 
278 
280 
298 


Int TRANSIENT rt 
ii TEMPERATURE ANDO POWER ACQUISITION PROGRAM FOR THE HEATER !11 
ttt AT THE BOTTOM OF CENTERAL COLUMN i 
OPTION BASE 1 

REAL Tul(35@) ,Yig(18@) ,Cas(350) ,Tem( 180) Vol( 350) ,Bath(3) Bat (3) 
CREATE BDAT "YIGIT: ,700,1" ,180 

ASSIGN @Path4 TO “YIGIT: ,70@,1" 

CREATE BDAT "“CANSU: ,700,1° ,350 

ASSIGN @PathS TO "“CANSU: ,700,1" 


Resist=1.886 
R2=Resist*Resist 
BEEP 


FOR I=l TO 350 
OUTPUT 709;"CONFMEAS OCV ,S@8-5@9,USE 0° 

ENTER 709;V,V1 

Vol(I)=¥ 

Tul(I)=sv!1 

NEXT I 

FOR I=! TO 180 

OUTPUT 709;"CONFMEAS DCU ,209,USE @” 

ENTER 709;T 

WAIT .09 

Tem(1)=T 

NEXT I 

FOR I=! TO 350 

IF I>180 THEN GOTO 308 

Yig( I )=.@006797+( 25825. 1328*Tem( I) )-(6@7789.2467#( Tem( I )*Tem(1)) )-(2195203 


4.3364¢( Tem( 1)°3) )+(8370810996. 1874*( Tem( 1 )*4))+.039 


300 
310 
328 
330 
348 
35 
368 
370 


Cas(Te(Vol( I )-Tul( 1) )/2.01e( Tul Cl (Vol CT Tul (1) 42.88 #0 4780/1147. )) 

NEXT I 

OUTPUT 709; "“CONFMEAS OCV ,317-319,USE 0° 

FOR I=t TO 3 

ENTER 709;Bat( 1) 

NEXT I 

FOR I=@ TO 3 

Bath(I )=.Q@006797+( 2582S. 1328*Bat( I) -(607789.2467*(Bat( I )#Bat(1)))-( 219520 


34.3364¢(Bat(1)°3))+(8370810996.1874*(Bat( I )*4))+.036 


380 
3398 
400 
418 
420 
430 
440 
45@ 
460 
470 
482 


NEXT I 

AS="TOP(C)” 
BS="MIO(C)* 
Cs="B0T(C)” 
PRINT USING “2X ,AAAAAAAAAAAAAAAAAARAAAA , 3X ,0D.00°sAS ,Bath( 1) 
PRINT USING "2X ,AAAAAAAAAAAARAAAAAAAAAAAR , 3X ,DD.00"3BS ,Bath( 2) 
PRINT USING “2X ,AARAAAAAAARAAAAAAAAAAAA , 3X ,0D.00°;C$ Bath 3) 
PRINTER IS CRT 

OUTPUT @PathS;Cas(«) 
OUTPUT @Path4;Yig(#) 
END 











APPENDIX D. TEMPERATURE ACQUISITION PROGRAM 


CCPECUP UEC ECE CCPECC EEC EEPOPEEPEG EEE PO CECE EOEC TE PerereeerCreeerertegyy 


1@ 
28 
38 
4a 
58 
68 
72 
8a 
98 
1@@ 
110 
120 
130 
140 
158 
160 
178 
180 
198 
220 
210 
220 
230 
240 


258 
260 
270 
280 
290 
300 


310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 


460 
470 
480 
4390 


tt 
1 


TEMPERATURE ACQUISITION PROGRAM vt 
FOR 15 HEATERS IN (STEADY STATE) ! 


CHEE ADL PS) CEE TE CEE EE TEE EEA 
OPTION BASE ! 

REAL Volts(18) 

REAL Temp(18) 

PRINT DATES< TIMEDATE) 

PRINT “INPUT AVERAGE POWER“ 

INPUT Power 

PRINTER IS 790! 

PRINT DATES( TIMEDATE ) 

PRINT “ 
PRINT Power , "WATT INPUT POWER“ 
PRINT " 15 HEATER ELEMENTS “* 


PRINT ” 3 MM SHROUQ SPACING 7 

PRINT “ “ 

PRINT “ ” 

PRINT * “ 

rerggeceecerercecegrerercnees SURFACE @1-5 
OUTPUT 709;"CONFMEAS DCV, 115-119,USE oe" 
FOR I#1 TOS 


ENTER 709;Volts(!) 
Tema I )=.0006797+( 25825. 1328eVolts(1)-(687789.24679(Volts( I evalts(1)) =< 
21952034 .3364e( Volts( 1 )°3) )+( 8370810996. 1874e(Volts(1)°4))+.039 


NEXT I 
reed 


Peeve reecuererreueel SURFACE #21-830 


OUTPUT 709; “CONFMEAS OCV 200-209 ,USE @” 

FOR I*6 TO 1S 

ENTER 709:Volts(!? 

Temp( 1 =.0006797+( 25825. 1328eVolts( I) )-(607789.2467¢(Volts( I eVolts(1)))-¢ 
21952034. 3364*(Volts(I 173) )+(8370810996. 1874*{ Volts(1)°4))+.033 


NEXT I 


tif PRINT OUTPUT TEMPERATURES 
PRINT * ° 
FOR I=! TO 1S 

PRINT “HEATER #"31,° TEMPERATURE IN DEG. C =";Temp(l) 


NEXT I 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


BATH TEMPERATURES (TOP TO BOTTOM)” 


OUTPUT 709;"CONFMEAS DCV 317-319 ,USE @° 

FOR I=16 TO 18 

ENTER 709;Volts(1) 

Temp( I )=.0006797+( 25825. 132BeVolts( 1) )-( 607789. 246 7#(Volts( 1 )eVolts(1)))-< 
21952034. 3364«(Volts( 173) )4+(83708 10996. 1874e(Volts(1)°4))+.036 

PRINT “TEMP. DEG. C “sTemp( 1) 


NEXT I 


PRINTER IS CRT 


ENO 
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APPENDIX E. POWER INPUT PROGRAM 


40 OUTPUT 709; "RST" 

Se@ OPTION BASE ! 

62 REAL Volts(!8) 

72 REAL Pow( 15) 

8@ REAL Ptot ,Pavg ,Perd:if 

‘90 OUTPUT 7@9;"CONFMEAS OCV 316 ,USE @* 

100 ENTER 709;Volts( 16) 

118 OUTPUT 7@9;"RST" 

120 OUTPUT 7@9;"CONFMEAS OCV 410-417 ,USE Q“ 

13 FOR I=! TO 8 

140 ENTER 709;Volts(I1) 

150 NEXT I 

160 OUTPUT 709;"“RST” 

178 OUTPUT 709;"CONFMEAS OCV,401 ,USE Q” 

180 ENTER 709;Volts(18) 

190 OUTPUT 7@9;"RST”" 

200 OUTPUT 709;"“CONFMEAS OCV 503-509 ,USE @° 

210 FOR J=9 TO {5S 

220 ENTER 709;Volts(1) 

230 NEXT I 

240 OUTPUT 709;"RST” 

250 OUTPUT 709; "CONFMEAS OCV ,S@@,USE 2” 

262 ENTER 729;Volts(17) 

270 OUTPUT 709; "RST" 

280 PRINT DATES( TIMEDATE ) 

290 PRINT * * 

300 PRINT ° * 

218 PRINT "INPUT VOLTAGE D.C. VOLTS PAN t";Volts(18) 

320 PRINT “INPUT VOLTAGE D.C. VOLTS PAN 2”°:Volts(17) 

330 PRINT ” “ 

340 PRINT ° * 

350 Resist#2.@ 

360 R2"Resist*Resist 

370 Pow( 3 )=(Volts( 18 )-Volts(3) )#(Volts(3)#(Resist+.55 )-Volts(18)¢.55)/R2 
382 Pow( 8 )=(Volts( 18 )-Volts(8))e(Volts(8)#(Resist+.S5 )-Volts(18)¢.S5)/R2 
390 Pow( 18 )=(Volts(17)-Volts( 10) #(Volts( 18 )*(Resist+.3)-Volts( 18)#.3)/R2 
400 Resist=1.972 

410 R2=Resist*Resist 

420 Resist=1.886 

430 R22Resist*Resist 

440 Pow ( tS #(Volts(14)-Volts(tS))/2.0te(Volts( 1S )-((Volts(14)-Volts(1S5))/2.01¢ 
4700/147.)) 

450 Resist=!.862 

460 R2=Resist*Resist 

472 Pow( 12 #(Volts(17)-Volts (12) e(Volts(12)e(Resist+.S)-Volts(17)*.5)/R2 
480 Powl 15 =(Volts(17)-Volts( 13) e(Volts(13).¢(Resist+@. )-Volts(17)#@. )/R2 
4908 Resist=!.852 : 

SQQ@ = R2=Resist*Resist 

S510 Pow(9)=(Volts(17)-Volts(9) )#(Volts(9)e(Resistt+.! )-Volts(17)#.1)/R2 
520 Resist#1.8182 

S30 R2Z=Resist*Resist 

S48 Pow(6 )=(Volts( 18 )-Volts(6))#(Volts(6 )e(Resist+.7)-Volts(17)°©.7)/R2 
$s Resist=1.8 

S60 R2=Resist*Resist 

S78 Pow(4)=#(Volts( 18 )-Volts(4) )#(Volts(4)*(Resist+.2)-Volts(18)*.2)/R2 
582 Pow(S )=(Volts(18)-Volts(S))#(Volts(S)e(Resist+.7)-Volts(18)*.7)/R2 
S98 Pow (14 )=(Pow( 1S )+Pow( 13))/2. 

FOO = Pout 7 )=( Pow( 6 )+Pow( 8) 72 
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Resist=!.6667 

R2=Resist*Resist 

Powll )=(Volts (18 )-Volts(1)e(Volts( 1 )#(Resist+.4)-Volts(18)e.4)/R2 
Pow( 2 )=(Volts(18)-Volts(2) )#(Volts(2 )e(Resist+.4)-Volts(18)#.4)/R2 
Powl1t)=(Volts(17)-Volts(11))e(Volts( 1} d#(Resist+.2)-Volts(17)#.2)/R2 
FOR I=! TO 15 

PRINT “BLOCK # ";I1," POWER = “;Pow(I)," WATTS" 

WAIT .S 

NEXT I 

Ptot=0 

FOR I=! TO 15 

IF I=t4 THEN 742 

Ptot=Ptot+Pow( 1) 

NEXT I 

Pavg=Ptot/14 

PRINT * 

PRINT "AVERAGE POWER = °;Pavg," WATTS” 

PRINT “ “ 

FOR I=! TO 15 

Perd: f#ABS( ( Pavg-Pow( 1) )/Pavg #100 

PRINT “HEATER #";:1," PERCENT DEVIATION “;Perdif 

WAIT .5S 

NEXT I 

END 
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APPENDIX F. TEMPERATURE ACQUISITION AND REDUCTION 


1Q 
20 
30 
49 
5Q 
60 
70 
80 
3a 
188 
110 
128 
130 
140 
1590 
160 
170 
18@ 
190 
200 
210 
co 
230 
248 
250 
260 
278 
52034. 
280 
2390 
302 
310 
320 
330 
348 
3550 
52034. 
360 
370 
380 
390 
420 


PROGRAM (STEADY STATE) 


PECHUCHE TERE C GE EEL EE EE ELE EOS EPEC CEEE EC CEP EE eee eed 
i TEMPERATURE ACQUISITION AND REDUCTION PROGRAM i 
Vt STEACY STATE (CENTER COL) tt 
't) PROGRAM PERFORMS THE REQUIRED FUNCTIONS NEEDED TO OBTAIN DATA it 
'! FROM THE ACQUISITION SYSTEM AND CONVERT THE THERMOCOUPLE QUT- tt 
rl PUT VOLTAGES INTO TEMPERATURE VALUES. i 
PUDUPP DATEL TEED ERP PEERS COPEL TERA T EE EEA EU PEE EPEC GPU T PTE 
OPTION BASE } 

REAL Velts(18) 

REAL T(18),Xtol 15) ,Nu(tS),Ra( lS) Delt (15) 

REAL Lnu( tS) ,bLgr( tS) 

REAL Power ,Tamb ,G,Kpg ,Oelx ,As ,Per ,Ta,Qconv ,Lbar 

REAL Qflux ,Tfilm,Beta ,Spvol ,Oynvisc ,Qcond 

REAL Kinvasc ,Kf ,Pr,Theta,H» ,Grm,Hux 

PRINTER IS CRT 

PRINT 

PRINT “ENTER AVERAGE POWER -WATTS” 

INPUT Power 

PRINT 

PRINT “ENTER SHROUD SPACING IN mm“ 

INPUT Space$ 

(PRINT * ° 

PEED UEC E ECE eee SURFACE #16-20 

OUTPUT 709; "“CONFMEAS DCV, 1t5-119,USE @° 

FOR I=1 TO 5 

ENTER 709;Volts(1) 

TCL =. 0@@6797+( 25825. 1328*Volts(1) )-(607789.2467#(Volts(I )eVolts(1)))-(219 
3364e(Volts(1)°3))+(8370810996.18749(Volts(1I)*4))+.039 

Xtol( 1 )=(31.75-6.985-( 1-16 )#1.27)/31.75 

NEXT I 

LHVT PEPE P PEEP eee teed anes SURFACE #21-%32 

OUTPUT 709; "CONFMEAS DCV 200-207 ,USE 9° 

FOR I*6 TO t3 

IF [#14 THEN 378 

ENTER 789;Volts(I) 

TCI )=.@006797+4( 25825. 1228 *Volts(1) )-(607789.2467e(Volts([)eVolts(1)))-€219 
3364*(Volts(1I)°3) )+(8370810996 .18749( Volts(1)°4))+.039 
Xtol(1)=431.75-6.985-( 1-16 )01.27)/31.75 

NEXT I 

OUTPUT 709;"CONFMEAS OCV ,209,USE @” 

ENTER 7@3;Volts( 1S) 

T( 1S )=,Q006797+( 25825. 1328eVolts( 15) )-( 607789. 2467#( Volts( 15 )*Volts(15)))- 


(21952034 .33642(Volts(15)°3) )+( 8370810996 .1874e(Volts(15)°4) )+.890 


410 
420 
430 
440 
52034. 
450 
460 
478 
460 


OUTPUT 709; "CONFMEAS DCV ,317-319,USE @° 

FOR I=16 TO 18 

ENTER 7@9;Volts(I) 

TCT =. @006797+( 25825. 1328eVolts( 1) )-(607789.2467#( Volts(1)*Volts(I)))-(219 
3364*(Volts( 19°35) )+( 8378810996. 1874*(Volts( 1 )°4))+.036 

PRINT “TEMP. DEG. C “;TCI) 

NEXT I 

PRINT “INPUT AMBIENT TEMPERATURE“ 

Tamb=(7(16 +T(17)4+T(18) 0/3 
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490 
500 
518 
520 
530 
548 
558 
560 
578 
S60 
598 
52 
618 
620 
630 
640 
650 
652 
670 
688 
6390 
708 
712 
720 
730 
740 
750 
762 
770 
780 
790 
B22 
818 
620 
830 
B42 
850 
860 
872 
880 
890 
900 
312 
920 


PRINTER IS 7@1 

PRINT “ * 

PRINT DATES( TIMEDATE ) 

PRINT “DATA REDUCTION OUTPUT" 


PRINT 

PRINT "POWER LEVEL (WATTS) :" Power 
PRINT 

PRINT "SHROUD POSITION (mm) 7" Spaces. 7 
PRINT 

PRINT “AMBIENT TEMPERATURE (C) :*,Tamb 
PRINT 

PRINT 

| 

t 

G=9.81 

Kpg=.1421 

Delx=.006731 


As=.0239* .0078 
Per=2+ .0239+2+.Q078 
Lbar=As/Per 
AS="Qcond" 
BS="X/L" 
D$="LOG(Nu- >” 
ES="LOG(Grxe)” 
FS="DELT” 
6S="H:." 
H$="Nu»" 
1$="Gréxe" 
JS="Temp” 
PRINTER IS CRT 


PRINT “INPUT CONDUCTION LOSSES AS COMPUTED FROM ELI.PACK* 

INPUT Qcond 

PRINTER IS 701 

PRINT “ CENTER COLUMN” 

PRINT " ” 

FOR Iel 70 15 

PRINT ” " 

PRINT " SURFACE &#" I 

PRINT " ° 

PRINT USING "2X ,AAAAA ,3X ,AAAR,3X ,AAAARAAA ,2X ,AAAAAAAAA 3X ,AAAA ,4X ,AA ,4X ,AAA 


1X ,AAAA ,2X ,AAAA” 5 AS ,BS DS ES FS 6S HS IS JS 


93@ 
7 940 
950 
960 
9370 
380 
9590 
1000 
1018 
1028 
1030 
1040 
1050 
1060 
1078 
1080 
1090 
1102 


Qconv=Power-Qcond 

Qflux=Qconv/As 

Delt (I )=T(1)-Tamb 

Tfilme273.154+(T(1 )+Tamb)/2 
Beta=(-7.9448SBE-Be(TFilm*2))+(5.7356479E-SeTfilm)-.0097610563 
Spvol=(4.6969SE-Se(Tfilm’2) )+(-2,.SZ74S5SE-GeTFilm)+.001341903 
Dynvisc#(3.23548511E-7e( Thilm’2) )4(-.00021474487*TFilm)+.036166792 
Kinvisc#SpvoleDynvisc 

Kfe(1.418181818E-3eTfilm)+. 1866 

Pr=(4.65706E~3*(Tfilm*2) )4(-2.922094*Tfilm)+465.3319 
Thetal=Powerl/( .O220«Kf ) 

Temp=T( I) 

Theta=0f lux*Lbar/Kf 

Hx@OQflux/Delt (I) 

Grm=(G*BetaeQfluxe(Lbar )°4)/(KfeKinvise*2) 

Nu» #Hx*Lbar/KF 

Nu( I )*Nux 

Ra( Il )#Pre6rm 
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11t@ 
1120 
113@ 
1140 
115@ 


Lor (I )=LGT(Grm) 

Lnu( I =LGT( Nux ? 

LgelLor (I) 

LynneLnu( 1) 

PRINT USING "1X ,0.0000,2x ,D0.00D ,2x ,D.DDD ,5X ,0.D0D ,4x ,D0.0D0D ,2x ,OD0D.0,1x,D 


6.000 ,00000.0 1X OU. .DD"sQcond ,Xtol( 1) inn bp Delt (1) Hx ,Nux ,6rm, Temp 


116@ 


117@ 
1180 
1190 
1200 
1210 
Vee 

1230 
1240 
1250 
1260 
1270 


NEXT I 


CREATE BOAT "NU: ,700,0,2" ,15 
ASSIGN @Path4 TO "NU: ,70@,0,2” 
OUTPUT @Path4;Nu(e) 

CREATE BDAT “RA: ,70€,0,2",15 
ASSIGN @PathS TO “RA: ,700,0,2” 
OUTFLT @PathS;Ra(s) 

PRINTER IS CRT 

CREATE BOAT “OTEMP2W3MM: ,720,17 35 
ASSIGN @Path? TO “DTEMPZWIMM: ,700,1° 
OUTPUT @Fath?;Delt(*) 

END 
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